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ABSTRACT 

The radiative acceleration of particles and the electrostatic potential helds that arise in low density 
plasmas hit by radiation produced by a transient, compact source are investigated. We calculate the 
dynamical evolution and asymptotic energy of the charged particles accelerated by the photons and the 
radiation-induced electric double layer in the full relativistic, Klein-Nishina regime. For fluxes in excess 
of 10 27 erg cm -2 s^ 1 , the radiative force on a diluted plasma [n < 10 11 cm~ 3 ) is so strong that electrons 
arc accelerated rapidly to relativistic speeds while ions lag behind owing to their larger inertia. The ions 
are later effectively accelerated by the strong radiation-induced double layer electric field up to Lorentz 
factors « 100, attainable in the case of negligible Compton drag. The asymptotic energies achieved by 
both ions and electrons are larger by a factor 2-4 with respect to what one could naively expect assuming 
that the electron-ion assembly is a rigidly coupled system. The regime we investigate may be relevant 
within the framework of giant flares from soft gamma-repeaters. 



Subject headings: acceleration of particles - radiation mechanisms: non-thermal 
gamma-rays: theory 



scattering 



1. INTRODUCTION 

Particle acceleration is an ubiquitous feature in many 
classes of astrophysical objects, such as active galactic nu- 
clei, gamma-ray bursts, pulsars. Despite several physi- 
cal processes are know to be responsible for accelerating 
particles under different conditions, radiative acceleration 
has special importance. Gurevich and Rumyantsev (1965) 
were the first to study the acceleration of charged particles 
caused by scattering off a radiation beam. They derived 
the relativistic equation of motion for the electrons and 
assumed that the radiative force is applied directly to the 
electron-ion assembly. For the typical radiative fluxes pro- 
duced in Type II supernovae, they found that ions may ac- 
quire energies of the order of their rest mass energy. Radia- 
tive acceleration was the subject of several further inves- 
tigations in the '80s, mainly connected with the possibil- 
ity that this mechanism could drive relativistic collimated 
outflows, or jets, such as those observed in Active Galactic 
Nuclei (AGNs) and in some Galactic sources (e.g. SS433). 
Earlier studies (e.g. Abramowicz and Piran 1980; O'Dcll 
1981a; O'Dell 1981b; Cheng and O'Dell 1981; Sikora and 
Wilson 1981; Phinney 1982) dealt with the radiative accel- 
eration of electrons assuming Thomson scattering in the 
particle frame and adopting a test-particle approach. This 
is fully justified when the medium is optically thin (for a 
discussion of the optically thick case see Turolla, Calvani 
and Nobili 1986) and allowed to clarify the role of Compton 



drag in establishing the limiting plasma velocity and the 
dynamical effects caused by anisotropic Compton losses of 
plasma internal energy ( "Compton rocket" ) . 

The first calculation of the radiative force acting on 
electrons taking into account the Klein-Nishina cross sec- 
tion was performed by Blumenthal (1974). More recently, 
Madau and Thompson (2000) re-addressed the issue of 
the dynamics of a cold plasma subject to an external ra- 
diation field in the Klein-Nishina regime. In a compan- 
ion paper (Thompson and Madau 2000), the feedback of 
pair creation on radiative acceleration was also considered. 
Contrary to earlier works, the main focus of these inves- 
tigations is on the dynamics of the fireball expected to 
energize gamma-ray bursts. In such a scenario, the mean 
photon energy is large enough to give raise to copious pair 
production. 

An interesting and related research thread is the acceler- 
ation by electrostatic fields generated by space charge sep- 
aration. Shvartsman (1970) considered a mechanism suit- 
able for positron acceleration in accreting neutron stars. 
He proposed that the radiation field incident on the plasma 
acts differently on electrons and ions because of the large 
difference in mass and cross section (cf. also Maraschi, 
Reina and Treves 1974; Maraschi, Reina and Treves 1978). 
Electrons are separated from ions and an electrostatic field 
arises to restore charge neutrality. The net electrostatic 
potential is able to accelerate positrons, which are limited 
by pair production and inverse Compton, but reach en- 
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ergies of tens of MeV. The calculation has been revised 
by Turolla et al. (1997), who found that the maximum en- 
ergy attainable is limited to <~ 400 kcV. Treves and Turolla 
(1999) and Treves, Turolla and Popov (2000) investigated 
this scenario in connection with the vacuum breakdown 
near an accreting black hole, while Hirotani and Okamoto 
(1998) studied the case of a supermassive black hole in an 
AGN. Litwin and Rosncr (2001) proposed another version, 
consisting of infalling blobs of plasma onto a neutron star, 
in which the electric field is generated by the co-rotation 
of the strong magnetic field of the star (cfr. also Hirotani 
and Shibata 2001). This mechanism has also several ap- 
plications to particle acceleration in solar flares and plane- 
tary bow shocks (cf. Ball and Melrose 2001 and references 
therein). 

In this Paper, we present a complete analysis of the ac- 
celeration of electrons and ions by electrostatic potential 
fields that arise in a low density plasma illuminated by a 
burst of radiation. Extending previous works (in partic- 
ular that by Madau and Thompson 2000), the dynamical 
evolution of the charged particles under the joint action of 
the internal electrostatic and the external radiation forces 
is investigated. In § 2 we derive the non-relativistic and 
relativistic equations of motion for the charged plasma par- 
ticles in presence of a Double Layer. Our results are pre- 
sented in §3 and thoroughly discussed in §4. Conclusions 
follow in §5. 

2. RADIATIVE ACCELERATION AND ELECTRIC 
RESTORING FORCE IN A PLASMA 

In a plasma, any slight deviation from a condition of 
global electric neutrality gives rise to large restoring forces. 
Radiation impinging on a plasma acts differently on ions 
and electrons, because of their large difference in mass, 
generating a space charge separation. Therefore, an elec- 
tric field rises to restore global neutrality. Since the re- 
sulting electric field is across two opposite space charge 
layers, it is often called "double layer" (DL), particularly 
with reference to studies in interplanetary bow shocks and 
solar flares (see e.g. Block 1978; Raadu 1989). 

Consider a transient, compact source of high energy 
photons of luminosity L surrounded by an optically thin 
shell of cold plasma. The plasma is composed of elec- 
trons and completely ionized ions of atomic number Z. 
The radiation flux is sufficiently strong that gravity can 
be neglected and, at the typical distances from the central 
source considered here, magnetic fields are assumed to be 
dynamically unimportant. The equations of motion of the 
electrons and ions (in a reference frame K at rest with re- 
spect to the central source) under the joint action of the 
radiative force and the restoring force of the DL are 

le-JJ (7e^eV e ) = Fe,c + F rad (1) 



dt 



' E,i 



(2) 



where 7 e ,i, ~v e ,i and m ej i are the Lorentz factor, 3- velocity 
and mass of electrons and ions, and F ra( j, F E , e and F E ,i 
are the spatial components of the radiative and electric 4- 
forces acting on electrons and ions. The radiative force act- 
ing on the ions is neglected because it is a factor (m e /rrii) 2 
smaller than that on the electrons. Here we consider mono- 
energetic particles. If the charged particles have a mo- 



mentum distribution, the previous equations must be av- 
eraged over the normalized phase-space density f e ,i(y e .i)- 
We checked a posteriori that the radiation-reaction force 
due to the energy radiated by the accelerated particles is 
negligible and does not influence the dynamics. 

The electric force is given by the spatial components of 
the Lorentz 4-force vector 

F E,e = -leZeE (3) 
F E ,i = JiZeE , (4) 
where E is the electric field and Z the atomic number of 
the ions. The general expression for the spatial part of the 
radiative 4-force acting on the electrons is (cfr. Blumen- 
thal 1974; Madau and Thompson 2000) 

'e , \ / v P • k \ I f da 



-H 



7e 



1 



dn 



-dedfldQ s 



(5) 

where e and k are the photon energy and unit propaga- 
tion vector, I e is the radiation specific intensity, da/dfl s 
is the differential scattering cross section, is the solid 
angle subtended by the source and fl s that over the 
scattering angle. The first term A(ek/c) = (eiki — 
ek)/c is the photon momentum change in a scattering 
(ei and ki are the photon energy and unit propagation 
vector after scattering), while the second 7 e [l — (v e • 
k/c](I € /e)(da/dil s )dedQ,dCl s represents the rate at which 
photons with energies between e and e + de and directions 
between il and f2 + d£l are Compton scattered by electrons 
with scattering angle between Q s and Q s + dQ s . The gen- 
eral expression for the differential scattering cross section 
is given by the Klein-Nishina formula 



da 
dSX 



16tt 




z i • 2 ' 
- - sin x s 



(6) 



where ot = 6.65 x 10 25 cm 2 is the Thomson cross sec- 
tion, primed quantities are evaluated in the reference frame 
K comoving with the electron and x s * s the (polar) an- 
gle between the incident and scattered photon directions. 
Equation (6) must be supplemented with the Compton 
formula for the energy of the scattered photon 

(7) 



1 + (e'/m e c 2 )(l — cosXg) 
The corresponding unprimed quantities in the system 
frame K are obtained Lorentz transforming photon en- 
ergies and directions by means of the the Dopplcr and 
aberration formulas. 

In the following we will investigate radiative acceleration 
from an isotropic, spherically symmetric source of radia- 
tion. Then, the motion of the charged particles is radial 
and the equations simplify to 
d 

le ~di ^ lemeVe > = E > e + Frad 



d 



(8) 
(9) 



1% ~dt ^ miVi ^ = Fe ^ ' 
where v e ^ are the radial components of the electron and 
ion velocities, and F ra d, F E , e and Fe,% are the radial com- 
ponents of the radiative and electric 4-forces acting on 
electrons and ions. They are given by 

F E , e = -leZeE (10) 
F E}i =nZeE (11) 

F rad = - I ^l_lH le (x _ k^-dednd^ X 2) 
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where E is the radial component of the electric field, and 
/i (/ii) is the cosine of the angle between the radial and 
incident (scattered) photon directions. In terms of these 

angles and the polar x s an d azimuthal <f> s scattering angles 

in K', it is \l x = // cosxl + (1 - M ) 1/2 sin Xs cos 4> s - 

If there are no external magnetic fields and the plasma 
particles move along the radial direction, the radial com- 
ponent of the electric field E can be obtained from the 
Gauss theorem. For a (geometrically thin) DL, we have 

E = AirZenAx, (13) 

where n is the volume density of charged particles and 
Ax = r e — Ti <C r 6j i is the displacement between the elec- 
tron and ion positions r e and r,. We checked a posteriori 
that the condition Ax <C r e ^ is satisfied. 

2.1. The non-relativistic equations in the Thomson 

regime 

The Newtonian equations for the radial motion of the 
electrons and ions under the joint action of the radiative 
force and the restoring force of the DL can be obtained 
from equations (8)-(12) neglecting terms 0(v/c). They 
read 

d 2 x e 



m e 



dt 2 

d 2 Xi 



= -ZeE + F, 



rad,T 



ZeE. 



(14) 
(15) 



where the displacements of electrons and ions from 

their initial position R (aj e ,i = fe,i — R)- For Thomson 
scattering (e 1 ~ e ) equation (12) gives 

F? a r d , T = c-T^ (16) 

where F = J I e /j,dedn = L/Anr 2 is the radiative flux emit- 
ted by the central source at distance r e = R + x e . 

By substituting equations (13) and (16) into equa- 
tions (14) and (15), we obtain 



d 2 3 



dt 2 m e 
d 2 Xi AwZ 2 e 2 n 



AirZ 2 e 2 n A a T F 
Ax + ■ 



m e c 



dt 2 



Ax . 



rrn 



(17) 
(18) 



These are the equations of motion of a particle oscillating 
with the plasma frequency (pulsation) 



' Ai{Z 2 e 2r r 



In particular, the electron plasma frequency is 

,1/2 , 



u P = 5.6 x 10 4 Z 



( n v 

VI cm- 3 / 



(19) 



(20) 



In a fully ionized hydrogen plasma (electrons and protons 
only), Z = 1 and rrii = m p . Finally, we can write two 
coupled equations for the displacements x e and Xi in the 
form 

-Ax+ ; C — n (21) 



(l + x e /R) 2 



Xi = ^-Ax , 
rrii 



(22) 



where 'x e ^ = d 2 x e> i/ d{ui e t) 2 and 
o~tL 



62 Z 



4:irR 2 cm e td 2 
n 



(r 



cm 



R 



10 9 cm 



. , cm (23) 
10 38 ergs -1 / 

represents the amplitude of the radiation-induced oscilla- 
tions (for small displacements x e <C R; see eq. [21]). 

2.2. The relativistic, Klein-Nishina regime 

The coupled radial equations for the motion of electrons 
and ions in the reference frame K are (eqs. [8] and [9]) 



7e^ ("f e m e Ve) 

(jirriiVi) = jiZeE. 



-JeZeE + Frad 



(24) 
(25) 



The radial component of the radiative force is calcu- 
lated inserting the Klein-Nishina and Compton formulas 
(eqs. [6] and [7]) into equation (12) and performing the 
integral over Q s . We obtain (see e.g. Blumenthal 1974; 
Madau and Thompson 2000) 



p — ° T 

r rad — 7< 
C 



/ 



K(x ) 



{ 7e 2 [Ml+/?e)-A(l+M 2 )] 



+ 



where (3 e — v e /c, x 
1 



x fi(l — jU/3 e ) j I e dedfl , 

7 e (e/m e c 2 )(l - I3 e jj) and 



(26) 



K(x) = 



0~T 

3 



2x 



+ 



Ax' 
-10a; 



2x' 



■ln(l + 2a; ) 



51a;' 3 + 93a;' 2 



51a; +9 



3(1 + 2a;') 3 



(27) 



(28) 



For Thomson scattering, equation (26) reduces to 

Frad.T - ^ [(1 + Pl)F - V e (U + P)] , 

where U — c _1 J I e dedfl and P = cr 1 J I e ^ 2 dedQ are the 
radiation energy density and pressure. The same expres- 
sion can be obtained Lorentz transforming the radiative 
force from the reference frame comoving with the electrons 
K . In that frame it is F rad — gtF jc and F is the radia- 
tive flux measured in the comoving frame. If the radiation 
field is stationary, the time component of the radiative 4- 
force is zero. Thus, in K it is F ra d — leF rad — "fe^rF /c, 
where F is given by the Lorentz transformation for the 
radiative flux 

F' ^^[(I + PDF-v^U + P)} , (29) 

and one finally obtains equation (28). 

For a point-like source or sufficiently far from a finite- 
size source the radiation field is in radial streaming, i.e. 
I e = J e ,o 6(n - 1) and F = cU = cP (in both K and K'). 
In this limit equation (28) reduces to 



a T F 3 



rad.T 



(30) 
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If, in addition, the input radiation spectrum is monochro- 
matic, I £j o = Io5{ € — eo), then equation (26) becomes 1 

~K(x' ) 



rad 



7 e TO e C 2 



(31) 



wherex = -f e (e /m e c 2 )(l— [3 e ). The term in square brack- 
ets represents the Klein-Nishina correction to the radia- 
tive force. Finally, substituting equation (31) into equa- 
tion (24) and using the expression for the electric field 
(cq. [13]), the equations of motion take the form 

(j e m e v e ) — —^irZ 2 e 2 nAx 



dt 



K(x ) 



1 + 



— (7 i m i t; i ) = 4TrZ 2 e 2 nAx . 
at 

Introducing the new variables y e 
tions (32) and (33) can be written 



-Ax + 7e 2 (l 



c 



(32) 

(33) 
equa- 



(1 



K(x ) 



1 + 



f x e /R) 2 



3 • m e A 
liVi = Ax , 

rrii 

where y e>i = dy e j/d(uj e t). 

3. RESULTS 



(34) 
(35) 



Assuming that x e ^ <C R and F = const (i.e. the lumi- 
nosity of the central source does not vary with time), the 
non-relativistic equations (eq. [21] and [22]) can be solved 
analytically. Taking as initial conditions x e ^(0) — and 
dx et i/dt\o — 0, the displacements of electrons and protons 
are given by 



x e (t) 



Xi(t*) X e 

where 

6 = s'm(uj r t - 



(1 + 6) 



4irZ 2 e 2 n 
<jtF 



u 2 t 2 - 



tt/2) 

m e + rrii 



(1 + 6) 



c 



m e mi 

1 + ^ 

rrii 



1/2 



1 + 



rrii 



(36) 

(37) 

(38) 
1/2 
(39) 

(40) 



cm e uj; 

The relativistic equations (eqs. [34] and [35]) have been 
integrated numerically by using a 4 t/l -order Runge-Kutta 
method (see e.g. Press et al. 1993). The initial condi- 
tions are the same as in the non-relativistic case. 6nce 
a solution for x e ^ is computed, the electric field of the 
DL is calculated from equation (13). We performed the 
integration for several values of the parameter C and for 
flashes with two different durations (r = 15.9(27r/u; e ) s 
and t = 4.8 x 10 4 (27r/cj e ) s). The results of our calcula- 
tions are shown in Figs. 1-7 for R = 10 9 cm, n = 1 cm~ 3 , 
Z=l. 

1 Equation (31) corrects a typo contained in a similar equation (eq.[48]) derived by Madau and Thompson (2000). 



4. DISCUSSION 

Equations (36) and (37) show the properties of the solu- 
tion in the non-relativistic regime. The system is a coupled 
oscillator with a driving force (radiation). When the ra- 
diation flash hits the plasma shell, electrons are pushed 
outwards and begin to oscillate. Because of the restoring 
force of the DL, ions are dragged by the electrons. The 
whole system oscillates essentially at the electron plasma 
frequency (w r /27r ~ o> e /27r) and is accelerated outwards. 
The amplitude of the oscillations of the ions is a factor 
m e /rrii smaller than that of the electrons. From equa- 
tion (36), it is easy to see that the mean acceleration 
imparted to the whole system by the radiation force is 
(ui 2 — co 2 )x eq = uj 2 x eq = <r es F/(m e + rrii), as one could 
expect. The amplitude of the oscillation of the electrons 
relative to the ions is < Ax >~ x eq ~ C Therefore, the 
average electric field of the DL is < E >~ AnZenC. In the 
non-relativistic limit, the velocity of the charged particles 
and the DL electric field increase linearly with C. 

Fig. 1 shows the numerical solution of the relativis- 
tic equations (eqs. [34] and [35]) for a low luminosity 
(L = 1.6 x 10 36 erg s _1 ), short duration (r = 1.8 ms) 
flash. The energy of the input monochromatic radiation 
spectrum is eo/m e c 2 = 0.01, sufficiently low that the 
Klein-Nishina correction to the scattering force is negli- 
gible. During the early evolution, the essential properties 
of the non-relativistic solution are recognizable. However, 
after a few plasma time-scales, oscillations are significantly 
damped. In fact, in the relativistic equations (e.g. eqs. [32] 
and [33] ) the radiative term contains a dependence on the 
electron velocity v e that gives rise to a force formally sim- 
ilar to a "viscous" force. However, this is not a physi- 
cal drag exerted by the radiation field onto the electrons 
because, by assumption, photons are streaming radially. 
This effect is simply a consequence of the Lorentz trans- 
formation for the radiation moments. In fact, the radiative 
flux measured in the comoving frame varies with respect 
to that measured in the system frame according to equa- 
tion (29). In radial streaming, the term v e (U + P) be- 
comes 2/3 e F. This accounts for the decrease (or increase) 
in the flux measured in the comoving frame caused by 
the electron motion and, ultimately, originates from the 
finite velocity of propagation of photons. This effect is 
already important at fairly modest electron speeds and, 
when v e > 0.01 — 0.1c, the decrease of the radiative flux 
makes radiative acceleration less efficient. The change in 
the sign of v e induced by the plasma oscillations causes a 
decrease/increase of the radiative force (as seen in the co- 
moving frame) that produces the same effect of a viscous 
force, damping the oscillations of the electrons on a few 
characteristic plasma time-scales (see Fig. 1) 

For C > 10 6 cm, the system approaches the truly rela- 
tivistic regime. As shown in Fig. 2, during the initial tran- 
sient phase, the radiative force is so strong that electrons 
are accelerated rapidly to relativistic speeds while the ions 
lag behind owing to their larger inertia. This produces an 
electric restoring force that accelerates the ions for all the 
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duration of the burst. Oscillations are completely damped. 
For a short duration flash, the final speeds and energies 
of the charged particles are rather modest. However, as 
shown in Fig. 3, for a long duration burst (r = 5.4 s), the 
asymptotic speeds approach the velocity of light. This is 
more clearly shown in Figs. 4 and 5, where the initial tran- 
sient is so rapid that v e ~ c in less than a characteristic 
plasma time-scale and the ions are later effectively accel- 
erated by the strong radiation-induced DL electric field up 
to 7 <~ 500. When the ions catch up with the electrons at 
t ~ 2000(27r/w e ), the electric restoring force becomes neg- 
ligible and the electrons are rapidly re-accelerated by the 
radiation field until their velocity becomes equal to that 
of the ions. In this asymptotic regime, the DL has already 
moved significantly from the initial position (r > 10 10 cm) 
and the Lorentz factor has become so large that the radia- 
tive flux in the comoving frame becomes negligible. Radia- 
tive acceleration is no longer very efficient and electrons 
and ions are almost freely streaming, with the electrons 
oscillating around the ions (see Fig. 5). The asymptotic 
large amplitude oscillations achieved by the electrons may 
give rise to strong internal dissipation, thereby increasing 
the plasma temperature. 

Models with C > 10 11 cm suffer, however, from two 
major limitations. In fact, assuming R ~ 10 9 cm, the 
luminosity exceeds 10 47 erg s _1 and the effective temper- 
ature of the radiation is > 10 9 K for a typical source size 
R s w 10 7 cm. The mean photon energy is then sufficiently 
large that the Klcin-Nishina correction to the radiative 
force (see eq.[31]) starts to become significant. This is a 
fortiori true if the primary photon spectrum has a signif- 
icant high-energy component. An example of the parti- 
cles dynamics in this regime is illustrated in Fig. 6 for 
monochromatic input photons with e /m e c 2 = 10. As 
expected, the main effect of the reduction in the effective 
cross-section is to lower the terminal Lorentz factor of both 
ions and electrons, but this does not change the essential 
features of the acceleration process discussed above. The 
second issue concerns the possibility that the wind may 
become pair-loaded as a consequence of 77 interactions 
among primary photons or between primary and once- 
scattered photons. This effect is bound to increase dra- 
matically the efficiency of radiative acceleration because 
it reduces the inertia of the plasma, as discussed in de- 
tail by Thompson and Madau (2000). Charge-separation 
in a pair-dominated medium is outside the scope of the 
present investigation, so in the following we restrict to 
cases with photon energies well below the threshold for 
pair production. For thermal radiation with effective tem- 
perature <5x 10 s K, the fraction of photons with energy 
two times larger than the threshold energy (~ 7 MeV, 
at which the cross section for pair-production peaks) is 
/ ss 10~ 54 . The luminosity from a source of characteris- 
tic size ~ 10 7 cm is L < 5 x 10 45 erg s _1 , corresponding 
to C < 3 x 10 9 cm for R = 10 9 cm. The optical depth 
for 77 interaction between two primary photons close to 
the source is t 77 rs 7i 77 ct 77 -R s , where cr 77 ~ <jt for large 
scattering angles and n 77 ~ (L/e)At/R^ for a burst of 
duration At. Although t 77 would be lower for beamed 
or anisotropic emission, close to an isotropic source it is 
much larger than unity and all photons above threshold 
can produce e+e~ pairs. However, for T < 5 x 10 8 K, 



there arc no photons with energy two times larger than 
the threshold energy (Ni w f(LAt/e) <C 1) and hence the 
radiation field is pair-free. As for the photons scattered in 
an off-radial direction by the plasma electrons at radius R, 
their optical depth to primary photons is ~ pi~ 77 , where 
p w nor Ax < 1 is the scattering probability. The value 
of pr 77 can be larger or smaller than unity, depending on 
the plasma density. But, even if pr 77 > 1, for T < 5 x 10 8 
K no additional pair-loading is induced because there are 
no primary photons above threshold and once-scattered 
photons did not gain appreciable energy in the scattering. 
Therefore, for C < 10 10 cm the dynamics of the plasma 
shells is not influenced by pair-loading. 

For these values of the parameter C, the probability that 
once-scattered photons collide with another electron is ac- 
tually larger than undergoing 77 interaction. If at the time 
of the second scattering event the electron has already ac- 
quired rclativistic velocity, in the electron rest frame the 
energy of the scattered photon is boosted by a factor 7. So, 
if p is the probability of the first scattering, the probability 
that a photon undergoes two scatterings is p 2 /7 (owing to 
the decrease in the Klein-Nishina cross section). Then, the 
number of twice scattered photons for a At — 1 s burst is 
Rj (LAt/e)p 2 /j ~ 10 5 V (£/10 46 erg s" 1 ). In order not to 
influence significantly the dynamics of the plasma shells, 
this number must be smaller than rts nR 2 Ax, which im- 
plies n < 10 11 cm~ 3 for Ax rs 10 7 cm. For plasma den- 
sities in excess of n ~ 10 11 cm" 3 , the dynamics of the 
plasma shells starts to be significantly influenced by Comp- 
ton up-scattering of off-radial photons (Compton drag). 
Twice-scattered photons are boosted to very large ener- 
gies. Energy conservation limits their maximum energy to 
7m e c 2 < 50 MeV, corresponding to a the total luminosity 
rj 10 4 V erg s" 1 . 

Compton drag caused by off-radial primary photons 
may also play an important role in the vicinity of the 
source. In order to understand qualitatively this effect, 
we solved the equations assuming a slight deviation of 
the radiation field from the perfect radial case (of the 
order of R 2 a /R 2 ~ 10~ 4 ). The solution is shown in Fig- 
ure 7 where, in order to avoid the regime of pair-loading, 
C ~ 3 x 10 9 cm. Compton drag decreases the electron 
Lorentz factor attainable during the electrostatic ion ac- 
celeration stage. Therefore, once the ions have been accel- 
erated to relativistic energies by the radiation-induced DL 
electric field, they overtake the electrons and the DL re- 
verses. At this stage the electrons are re-accelerated at the 
expense of the electrostatic potential energy of the DL and, 
ultimately, of the kinetic energy of the ions. The asymp- 
totic Lorentz factors of ions and electrons are about 10 
times smaller than those achieved for a perfectly radial ra- 
diation field (point source). Placing the plasma shell closer 
to the central source would end up in a sufficient deviation 
from radial streaming to cause too a large decrease of the 
overall efficiency of the acceleration mechanism. 

We note that, even if the plasma is initially cold, random 
internal energy will be imparted to the electrons by the in- 
cident hard-X and gamma-rays, because direct energy and 
momentum transfer results when outward-directed pho- 
tons are removed from nearly radially streaming radia- 
tion field and scattered away. In principle, if a sufficient 
amount of internal energy is transferred to the electron 
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plasma, it can be effectively converted into bulk motion 
through the so-called "Compton rocket" (O'Dell 1981a; 
O'Dell 1981b; Cheng and O'Dell 1981). In this way, the 
energy lost to heat may be compensated by that gained 
by the electrons thanks to the additional Compton thrust. 
Although the study of this regime is beyond the scope 
of the present investigation, the outcome of the electron 
acceleration process with the additional Compton thrust 
may not be qualitatively different with respect to what 
found here. 

The radiation-induced electrostatic acceleration of the 
plasma shells is effective if they are quite spatially con- 
centrated and overdense with respect to the surrounding 
medium. In fact, if the shells are spatially extended, the 
escaping radiation would accelerate the plasma located at 
larger and larger radii at progressively lower speeds. The 
innermost, high velocity DLs would rapidly overtake the 
outer, low velocity ones, producing shocks and dissipation 
of kinetic energy into heat. 

It is noteworthy that, for intense radiation fluxes (C > 
10 8 — 10 9 cm), the asymptotic energies achieved by both 
ions and electrons are larger by a factor 2-4 (depending 
on the burst duration, see Fig. 7) with respect to what 
one could naively expect assuming that the electron-ion 
system was rigidly coupled with a total mass (m e + m;) 
and electron scattering cross section gt (see e.g. Gure- 
vich and Rumyantsev 1965). This is caused by the fact 
that the radiative force is not conservative and hence 
the dynamical evolution of the system is of fundamen- 
tal importance in determining the asymptotic energies. 
Thanks to the prompt response of the electrons, the work 
done by the radiation field is accumulated in the DL and 
is later effectively converted into kinetic energy of the 
ions through the acceleration produced by the DL elec- 
tric field. As shown in Fig. 7, this process is significantly 
more effective than radiatively accelerating a rigidly cou- 
pled electron-ion system. In fact, ions become relativis- 
tic at t\ ~ 30(27r/cj e ) and are electrostatically acceler- 
ated to large Lorentz factors up to i 2 — 300(27r/cj e ), 
when they catch up with the electrons. In the interval 
t\ < t < t 2 , the average velocity of the electrons is smaller 
then that of a rigidly coupled system and hence the ra- 
diative force is more effective. The net work done by the 
radiation field is 2-4 times larger than it would otherwise 
be and is stored in the DL as electrostatic potential en- 
ergy. The maximum electric field and energy can be es- 
timated assuming that, initially (0 < t < t\), the elec- 
trons travel at the speed of light while the motion of the 
ions is approximately non-rclativistic. Then, x e ~ ct and 
d 2 Xi/ d{oj e t) 2 ~ (m e /mi)x e = (m e /m,i)ct. Integrating this 
equation with Xi(0) — and dxi/dt\ = 0, we can calcu- 
late the time t\ at which Vi ~ c and the separation Ax max 
at t\. It is: ti = (2m i /m e ) 1 / 2 (l/uj e ) ~ 10(27r/w e ) and 
Ax max — (2\/2/3)c/wi — 2 x 10 7 cm. Substituting these 
expressions in equation (13), we obtain 



En 



8ttV2 Zenc 



/ n \ 1/2 

0.13 1 -) stVcm 

VI cm - ' 3 / 



(41) 

in quite good agreement with the maximum electric field 



reported in Fig. 7. The distance traveled by the ions dur- 
ing the relativistic acceleration phase {t\ < t < t 2 ) is 
~ c(t 2 — h) and then their asymptotic energy is approxi- 
mately given by W max ~ eE max c(t 2 -h)~ 10 11 cV. 

Despite the large asymptotic energies achieved by ions 
and electrons, the overall efficiency of this process is very 
low. The final energy acquired by an ion is E tot ~ 
ZmpC 2 -foo ~5x 10~ 2 Z (700/100) erg. Comparing this en- 
ergy with the total work done by the radiation field on Z 
electrons, W rad ~ ZLtct t /AttR 2 - 5 x 10 2 Z (L/10 46 erg) 
erg, we find that the efficiency is not larger than <~ 
10- 4 ( 7 co/100). 

5. CONCLUSIONS 

We have shown that radiation-induced electric fields can 
increase significantly the terminal energy of a ion-electron 
plasma accelerated by a transient radiation source. The 
asymptotic Lorentz factors are a factor 2-4 larger than 
those attainable by radiative acceleration alone. The ef- 
fect is maximum when both radiation drag and 77 interac- 
tions arc negligible, i.e. when the central source is distant 
(R 2 /R 2 < 10~ 4 ) or beamed, the plasma density is low 
(n < 10 11 cm~ 3 ) and the average photon energy is well 
below the threshold for pair production (T < 5 x 10 8 K). 
For R ~ 10 9 cm, these limitations imply L < 5 x 10 45 
erg s _1 , or C < 3 x 10 9 cm. Decreasing/increasing the 
initial position of the plasma shell by an order of mag- 
nitude results in a corresponding decrement/increment of 
the maximum allowed luminosity (L < 5 x 10 43 erg s _1 
and L < 5 x 10 47 erg s _1 , respectively). However, because 
C remains the same and, for a fixed duration burst, is 
proportional to the total work done by the radiation field 
on the ion-electron assembly, the dynamics and terminal 
energy are similar in all these cases. 

Although direct astrophysical applications are outside 
the scope of the present investigation, it is quite interesting 
to note that giant flares from soft gamma repeaters (SGRs) 
are typically within the parameter range considered here. 
The total energy emitted in two giant flares observed from 
SCR 0525-66 and SGR 1900+14 was w 10 44 - 10 45 erg on 
a timescale w 1 s, with an implied luminosity of about the 
same order (see e.g. Hurley 2000 for a review). In the 
picture proposed by Thompson and Duncan (1995) the 
giant flare is produced by a large-scale instability of the 
neutron star magnetic field which leads to the formation 
of an expanding pair fireball. The radiative acceleration 
process presented in this paper has no direct connection 
with the dynamics of the baryon-shcathed pair fireball, 
whose characteristic size is w 10 6 — 10 7 cm. However, it 
is possible that the outflowing radiation front impinges on 
a low-density ordinary plasma at much larger distances, 
w 10 9 cm, accelerating protons and electrons to relativis- 
tic velocities as illustrated e.g. in Fig. 7. 

This work has been partially supported by the Italian 
Ministry for Education, University and Research (MIUR) 
under grants COFIN-2000-MM02C71842 and COFIN- 
2002-027145. We thank an anonymous referee for helpful 
comments. 
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Fig. 1. — Displacement x, velocity /3 and Lorentz factor 7 of electrons (solid line) and ions (dashed line) for £ = 10 5 cm and r = I5.9(2iz/ui e ) 
s. For R = 10 9 cm, n = 1 cm -3 and Z = 1 this corresponds to L = 1.6 X 10 36 erg s _1 and r = 1.8 ms. The induced DL electric field E is 
also shown. Here the energy of the monochromatic radiation spectrum (eo/m e c 2 = 0.01) is sufficiently low that the Klein-Nishina correction 
to the scattering force is negligible. 
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Fig. 4. — Same as Fig. 1 for C = 6 X 10 s cm and r = 4.8 X 10 4 (27r/aj e ) s. For R = 10 9 cm, n = 1 cm -3 and Z = 1 this corresponds to 
L = 9.7 X 10 44 erg s — 1 and r = 5.4 s. The scale on the x-axis is logarithmic. 
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Fig. 5. — Same as Fig. 1 for C = 10 12 cm and r = 4.8 X 10 4 (2x/iO e ) s. For _R = 10 9 cm, n = 1 cm" 3 and Z = 1 this corresponds to 
L = 1.6 X 10 48 erg s — 1 and r = 5.4 s. The scale on the x-axis is logarithmic. The dotted line represents the motion of a single particle with 
total mass (m e + mi) and electron scattering cross section <tt- 
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Fig. 6. — Same as Fig. 1 for £ = 10 12 cm and r = 4.8 X 10 4 (2x/iO e ) s. For _R = 10 9 cm, n = 1 cm" 3 and Z = 1 this corresponds to 
L = 1.6 X 10 48 erg s — 1 and r = 5.4 s. Here the energy of the monochromatic radiation spectrum is assumed to be €o/m e c 2 = 10. The scale 
on the x-axis is logarithmic. The dotted line represents the motion of a single particle with total mass (m e + m,;). 
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Fig. 7. — Same as Fig. 6 for C = 2.8 X 10 9 cm and r = 4.8 X 10 4 (27r/w e ) s. For R = 10 9 cm, n = 1 cm -3 and Z = 1 this corresponds to 
L = 4.5 X 10 45 erg s -1 and r = 5.4 s. Here the anisotropy of the radiation field has been accounted for (see text for details). 



